EFFECTS OF SUCTION AND INJECTION ON SELF-SIMILAR 
SOLUTIONS OF SECOND-ORDER BOUNDARY-LAYER EQUATIONS 



By 

SHABBIR AHMAD 


Me 

\Sl3 

M 




TH; ■ , 

n(5 


PiK s 1,^ 



DEPARTMENT OF MECHANICAL ENGINEERING 


INDIAN INSTITUTE OF TECHNOLOGY KANPUR 





EFFECTS OF SUCTION AND INJECTION ON SELF-SIMILAR 
SOLUTIONS OF SECOND-ORDER BOUNDARY-LAYER EQUATIONS 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 



By 

SHABBIR AHMAD 


I'H tT 



Mill 


to the 

DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

JULY 1973 


me- 





CERTIFICATE 



Certified that this work on "EFFECTS OP SUCTION AND 
INJECTION ON SECOND ORDER SEIF-SBIILAR SOLUTIONS OF 
BOUNDARY LAYER FLOWS" hy Shabhir Ahmad has been carried 
out under my supervision and that this has not been 
submitted elsewhere for a degree. 



Noor Afzal 


Assistant Professor 
Department of Aeronautical Engineering 
Indian Institute of Technology Kanpur 







■'^.^.73 2jf 





Ill 


ACKKQ'gLEDSKTvTffT^-'P.q 


I an highly indebted to Dr. Koor Afzal for 
a^gestingthe problem and for his guidance throughout the 
course of this v/ork. 

My thanks are also due to Mr. S.C. Raisiiighani 
for his help during the present work, 

I also thank Mr. G.l. Misra for typing the 
manuscript neatly and efficiently. 


SHiBBIR mUH 





SYNOPSIS 


Effect of sdction or injection on self-simlar solutions 
on second order (due to longitadinal cuirvature, transverse 
curvatare, displacement speed, external verticity and temper- 
ature gradient) boundary layer flows has been studied. Two 
cases of prescribed wall temperature and that of insulated 
wall with full similarity in viscous dissipation are 
considered. Numerical solutions have been obtained for 
large values of parameters and the results are critically 


discussed. 





CHiPTER 1 


IN'ERODTJCTIOR 

The problem of boundary layer flows with suction or injection 
is of great importance in engineering applications, like boundary 
layer control, transpiration cooling ard other diffusional operations. 
It has been studied extensively and references thereof can be found 
in the texts by Schlichting (I968), Eckert (I972) and Rosenhead (1965). 
These studies have employed the classical boundary layer equations 
and hence are valid for very large values of Reynolds numbers. It 
moderately large Reynolds number, however, it is wellknown that the 
boundary layer equations require certain second order corrections, 

_i 

(VanDyke I962) whose order of m^nitude isR ^ when compared to 
Prandtle's boundary layer. These corrections may be attributed to 
arise from longitudinal curvature, transverse curvature, displacement 
speed, external vorticity and stagnation enthalpy gradient. 

In the present woifc we study the effect of suction or 

injection on these second order corrections. To this end we have 

studied the effect of suction and injection on self-similar solutions 

of second order momentum and energy equations with full similarity 

in viscous dissipation. 

* 

For an impermeable wall tfce general structure of self 
similarity for second order effects have been studied by Afaal 
and Oberai ( 1972 ) and Warle and Davis ( 1976 )* These authors have 
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shown that in addition to Palkner-Skan pressure gradient parameter g , 
the second order equations contain three new parameters , due to 
longitudinal curvature, due to transverse curvature and due 
to displacement effect. Further it is shown that for a given g , 
the second order solutions show a large numher of singularities for 
certain critical values of A^, and A^, Physical as well as 
mathanatical structure of these singularities is explained by 
Raisin^ani ( 1973 )* particular it is shown that the singularities 
occur whenever the pammeters A^ etc. assume a value equal to eigen 
value of the c'orre spending homogeneous problem with homogeneous 
boundary conditions. 

Generalization to heat transfer with full similarity in 
viscous dissipation has been made by Raisinghani (1975) and Afzal 
__and Raisinghani (1975) • They have studied the tr/o cases of prescribed 
Wail temperature and of insulated wall. For the latter case the 
classical concept of recovery f aotor has been extended to second 
order flows. For longitudinal curvature problem, the jointly self- 
similar solutions for second order momentum equations have been 
studied by Harsimha and Ojha (1967)> while the corresponding heat 
transfer problem for a prescribed wall temperature has been 
studied by Kadambi and Gupta (1971)* 

For our study, the self -similar solutions with suction 
and injection have been obtained following the analysis of 
Raisinghani (l975)j which involve an additional suction parameter G. 
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(G > 0 represents the suction while C < 0, the injection). Bie 
equations have been integrated numerically for various values of C 
in the raa ge of -.6 -^,6. To see clearly the effect of suction 

and injection, solutions are obtained for g == 0, .5 and 1, (as they 
represent important cases of flat plate, plane stagnation and 
axi symmetric stagnation flows), Prandtl number c = .7, 3 and 5 > and 
for a large values of etc, in the mnge of "7 *5 £ , A^, Jl 4* 

The transverse curvature problem with suction and injection 
has been studied earlier by Wannous and SparTOw (1965)» These aathors 
have studied a particular problem of momentum transfer in axial flow 
along a circular cylinder and have shown that suction reduces the 
second easier skin friction while injection accentuates it. In tuamis 
of cur parameters their analysis oora:«sponds to 3 = 0 and A^ = 0. 






CHAPTER 2 


FIRST MD SBCOHI) ORDER EOulDARi LAYER EQTJAT'IOHS 

2 . 1 Introduction 

The second order houndary layer equations sxe obtained from 

Navier-Stokes equations using the method of matched asymptotic 

expansions, VanDyke (1962). In this method the flow field is divided 

1 

into two regions me near the wall of the order of R ® and the other 
outside this region. The two regions are described by two separate 

_^i 

asymptotic expansions with R ^ as an expansion parameter. The two 
expansions are matched in the overlap region. 

2 .2 Pormulation Of Second Order Problem 

The Navier-Stokes equations in a non-dimensional vector form 
(here distance X(s,n), velocity U(U,Y) are non diiientionalized wife 

respect to characteristic length L of the bo -.y and reference velocity 

2 2 

respectively, pressure by pD„ and temperature byU„/G^) maybe 

G V p 

written as 

div TJ = 0 , (2.1) 

U. grai U + grad P => -R curl curl U (2.2) 

U. grad T - r"'' T = r”^ grad U. def U, (2.3) 

where R = Uc Lp/v is the characteristic Reynolds number, 0 = V 

the Prandtle number. In writing the above equations, the flow is 
assumed incompressible. 




The houndary conditions at the surface are 


(2 4) 
(2.5) 


U = 0 , V = v^(s) 

T = T^(s) or grad T = 0 

where V^(s) is of the order c£ S' ‘ J? upstream the flow has to 

approach prescribed, may be non-uniform, velocity and temperature 
fields. 


The outer asymptotic expansions in the limit X(s,n) fixed 
R CO are 


B = TJ^ (s,n) + a ^ TJ 2 (s,n) + ... 

(2 .6a) 

1 

1 = T^(s,n) + R ® + ... 

(2.6b) 

P = P^(s,n) + R ® PgCsjn) + ... 

( 2 . 6 c) 

T = T^(s,n) + a ^ T 2 (s,n) + ... 

(2.6d) 


Substituting these expansions in (2, 1-2. 5) and collecting terms of 
equal powers of S we get the well-known Euler’s equations whose 
integration along stream line gives 

+ lJ^/2 = (2.7 a) 

T^ = (2.7b) 

where is the first order outer stream ftmotion and and 

H^(f^) are Bemoullis functions to be determined for upstream 


conditions 
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The secoid. order equations are pertarhation form of Euler 
equations and their integration gives 

Pg + Ui , hg = ^2 

Tg « ’i'2 

If the on-coming stream is independent of Reynolds nunher, then 

B 2 “ 0 = 

1 

The outer equations are invalid in a region of order R near 
the wall. To study flow near the wall, the appropriate (Prandtle) 
variable 

I = n R® (2.9) 

a nd the corresponding inner limit K fixed, R gives the following 

inner asymptotic expansions 


1 

U = u^(s,N) + R ^ UgCsjU) + 

(2,10a) 

1 -1 

¥ = R 2 v^(s,N) + R V2 (s,N) + 

(2.10b) 

P « p^(s,]j) + R ® PgCs,]?) + ... 

( 2 .IO 0 ) 

1 

T =. t^(s,N) + R ^ t^Cs,!!) + ... 

(2.10d) 


Using these expansions in equations (2. 1-2. 5), we get first and secmd 
order governing equations for inner (boundary l^er) flow, 

_i 

These inner equations are valid in a region of order R ® and 
violate the conditions at infinity. Matching the two solutions we 
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obtain the missing boundary conditions needed to solve the boundary 
layer equations. 

2.2.1 First Order Boundary Layer IBquations 
Gontinui ty 

(2^^ ^1)3 ° » (2.11) 


Momentum 


’'ls(=>0) -“iM - 0- 

Energy 


"1 *1H ’'1 *1H - - 0 • 


Boundary and Matching conditions are 
u^(s,0) = 0 , T^(s,0) =. \(s) 

u^(s,H) = U^(s,0) as E ^ « 
t^(s,0) = t^(s) or t^j^(s,0) = 0 

t^(s,N) = H^(0) as E . 


(2.12) 


(2.15) 


(2.14) 


( 2 . 15 ) 


2.2.2 Second Order Boundary Layer Equations 
Continuity 


^^^'^2^3 “ "• r ^^(0 00s e/r) E u^ ] 


' 2 ''E 


1 s 


- [ (K + 0 00 s e/r) E v^ 1 ^ , 


( 2 . 16 ) 
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Momentum 

^ 2 s ^2 ^ 1 s ^ 2 N "^2 ^ 1 ir " ^ 2 M 

= [ N - N U^(s, 0 ) U^g(s, 0 ) + 

- [KK (s, 0 ) + K r {IJ? (s,0) - }dN]g + 0 cos e/r 

- B^(0) VgCsjO) + I U^(s,0) UgCs.O)]^ (2.1?) 

Energy 

^ 1 ^ 2 s ^^2 ■^ 2 ^ 1 N ” ^ "* W " ^ 2 N 

= K (N t^g + t^jj - 2u^ u^jj) + 3 cos e t^jj/r (2.18) 

Boundary and matching conditions are 


u^CsjO) = 0 , T^CsjO) » 0 

UgCsjN) = H KU^(s, 0 ) + r^ N B’Cc) + U2(s,0) as S « 

tgCsjO) = 0 or tgjjCsjO) = 0 


( 2 . 19 ) 


(2.20) 


t2(s,lT) = H!j(0) ^^(s,ll) 1 “ 

3Hi 

^hereH!^( 0 ) •= (^}v {,_^,0 

The second 03 :?der equations are linear and hence can he 
divided into a number of simpler problems, Bae terms proportional 
to K, 3 cos e/r , TJ 2 (s, 0 ), B 1 >( 0 ) and Htj(o) may be attributed to 
arise due to longitudinal curvature, transverse curvature, displacement 
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speed, external vorticity and temperatare gradient respectively. 

In the following pages qpsn tit ies ( s,0) , U^CsjO), H^(0), 
(o) will he written without their arguments as and B^ . 





CHiPTSE 3 


.^^LP-SIMIllR iJjiLYSlS 


3.1 Introduction 


As mentioned eaxlier, self -similar solutions of first and 

second order 'boundary layer equations have 'been obtained “by Afzal 

and Raisinghani (1973) for flow over an impemea'ble surface (V =0) . 

w 

Here we shall study self similar solutions of the first and second 
order problem in the presence of suction and injection at the wall, 
for an incompressible two dimensional and ajdsymmetric flows. 


The governing equations for first and second order problan 
are transformed using the Gortlers variables 


fS 

= / H 

0 


r^^ ds 


r^ N H 


1 




( 3 . 1 ) 


The equations obtained in terms of ? and n axe subjected to the 
condition of similaxity and finally equations with full similarity 
in dissipation are obtained. 


3 *2 Governing Equations For Hon-Sindlax Flows 
3 . 2.1 Hirst Order Problem 


The first order stream function 1 })^ defined by equation of 
continuity ( 2 , 9 ) 






-r'' V, 


9 


(3.2) 
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The stream function and temperature t^ are expressed as 


tI)^(s,n) = ^ £C?,n) , 

t^(s,N) = (t^ - gC5,n3 + Hj. 

Using (3*2) and ( 5 . 3 ) we get for the velocity components 
= U^ f 

= -Uj r^[£ + + CAy + j A^-l) n£’]/v^ 

where A is defined hy 
r ^ 


A 


<!> 



d? 


The governing equations (2.12 - 2.15) may he written as 
momentum 


( 5 . 3 ) 

(3.4) 

(5.5) 

(3.6) 

(5.7) 


ftit + ££,. + A _ 25(f'f’ - f^£”) = 0, 


Energy equation 


a“^ g" + f g’ - f’ g - 2eCf' gg - £^ g’) = - 


E 


V 1 


Boixndary conditions 


(3.8) 


(3.9) 


f(5,0) + 2S £^CC,03 = c(0 


fCe,o) = 0 , = 1 


(3.10a) 

(3.1 Oh) 


gC§,0) = 1 , gC5,“3 = 0 . 


(3.11) 


Where 0(5) = - j^/(U^ r^) is the suction parameter and 

2 

E = U^/(t^ ~ 2^) is the Eckert mmher, suffix C and dash represents 

differentiation with respect to 5 and n respectively. A and 

1 

A, TT a-re defined hy (3*7) and are known as pijincipal velocity function 
w 1 
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principal thermal function respectively. These principal functions 
are of fundamental importance in describing flow conditions. In 
non-similar flows, the principal functions vary from point to point 
while they remain constant along a stream line in case of similar 
flows. Moreover if the two flows are similar, the value of 
principal functions remain the same for the corresponding stream 
lines. 

(5.2.2) Second Order Problem 

The second order stream function \|)2» satisfying equation 
of continuity (2.14), is defined as 


^ 2 N ~ ^^^"2 3 N u, cos e/r], 

(5.12) 


'^2s ~ 3 cos e/r) N v^j. 

In terms of variables C , h the stream function and temperature 
t^ are expressed as 


FC 5 ,n) 


(3.13) 


t2 = Ct,, - Hi) GCg,r,) (3.14) 

The velocity components obtained from (3. 12) and ( 3 . 13 ) are 

U 2 = Uj [F’( 5 ,n) - p £•] (3.35) 

V2 = r 3 [-{F + 25 F^ + + 3 - D F' } 

+ p(B^ + B^) {£ + 25 + pCAu + 3 A^ - 1)£‘}]//2C (3*16) 

where Bjj, and are defined by (5.22) and (3.25). Using 

(5.12-5.16) and the second order governing equations (2.15-2.18), '®'e have 
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Momentum equation 


F’" + £F" - 2Au £»F' + £"F - 2? F' “ f" F^+£»F^ - f^F”) 


BJ{-tiC 1+A„ ) £*” + (A + Ay - 1) Cf” + ff’ + 25 f' 

“z 1 


+ H r + 2Ay ) (n Ay + a + 2? a^) } 

n 5 £ 1 1 


/Cl + A )] + Bt[-TiC2A„ + £’”) + + £f* - Ag n 

1 1 t 


•pt2 


+ 25 £^ £'] - B^Ca + 25 a^) - B^C2Ay^ + Ag^D 


Energy equation 


a"^ G” + £ G' - A^ _y f’ G - 25 (£' G - £ G’) 

w~”l ^ 


^ _y . g - F g* - 2CCF’ - F^ g'} - 2E £” F' 


= A F’ 

t -H, 
w 1 


+ Cn g')' + E F' C-n f' + 2£')] + B^[-a“-^Cn g'] 


- 1 , 


+ E £" Cn F' + 2£’)] 


Botmdaxy conditions 

FC?,0) + 25 FgC5,0) = 0 
F>C5,0) = 0 

F’Ce,n) = n(-B^ + ^t V ®d F “ 
GC5,0) = 0 , G'C5,0) = 0 


(3.17) 


(5.18) 

(3.19a) 

(3.19b) 

(3.19c) 

(3.20a) 


G(5,n) = Bg(n-a) as n -»■ “ 
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where -various quantities are defined hy 

a = lim Cn-f), (5»21) 

^ 00 

K/Cr^ up, = -/n 3 cos e/Cr^J Up (5.22-2p) 

Bd = ByUj , 3^ = Bpu2 (3.24-25) 

Bg = ^ HJ/Ct^ - Hp, (5.26) 


, A-o , are principal functions of longitudinal 

^ t d 

curvature, transverse curvature, and displacement speed respectively 
and are defined hy (3.7)* 

3.3 Self Similar Solutions : 

3.3*1 First Order Problem 

To study self -similar solutions of first order boundary layer 


equations, we assume 


f = fCn) , 

(3.27) 

and A„ to be a constant say g . 

1 

(3»10) reduce to 

Now the equations (3.S) and 

f’" + £f" + 6Cl-f’^) = 0 

(3.28a) 

£(0) = C 


£' (0) = 0 , £‘ C“3 = 1 

(3.28b) 


where G, the suction parameter is to be constant. 

linearity of first order energy equation permits its study 
independently of Eckert number, i.e,, 

g(5,n3 = giC5,n) + E g2CS,n) 
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If and are similar, they must be functions of n only. 

But still E may be a function of C , which renders g non- 
similar. Hence for full similarity with viscous dissipation 

E has to be constant, liiich implies g. = 2 g , a constant, 

w~ ^1 

Therefore for full similarity, we assume 

gCn) = + E g2Cn) (5.29) 

substituting (5.29) into (5.9) an-<3. collecting terms proportional 


to powers of E, we get 

g'l + £ ~ 2B f« = 0 , (3.50a) 

g^(o) = 1. gi(») = 0 , (3.30b) 

g" + £ gl - 26 f' g- + f”2 = 0 , (5.31a) 

g2C0D=0=g2W. ( 5 . 51 b) 

3.3.2 Second Order Problem 


The second order boundary layer equations are also linear 
and therefore allow superposition. These hi^er order equations 
can be divided into a number of simpler problems. In the 
equations (5*17) and (5. 18), the terms proportional to arise 
due to longitudinal curvature, due to transverse curvatuiE, 

B^ due to displacement, B^ due to external verticity and due 
to temperature gradient in the on-coming stream. 

Further the second order heat transfer pix>blem like that of 
first order can also be solved independent of Eckert number. 
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Therefore to obtain the self-sisdlar equations for the second 
order problen we asaime 


I fW (n) . 

m=£,t,d,v 

(3.32) 

1 [gJ®^Cp) + EG2^®^Cn)] 

m=£,t,d,v,e 

(5.33) 


’ '^t ’ 

Using equations (3.32-5*53) a^d separating the various second order 
effects, the equations (3. 17 “3 *20) take the form (using operator 
notation) 



are constants say , A^ and A^ respectively. 


F»'i + f pn . (20 + xr) £' F‘ + (1 + X^) f”F = % (5.34a) 

FfO) = 0 = F’CO) , FC») = (3.34b) 


+ £G’ - (2&+Xi)£»G^ + Cl+H)Fg^ - 23 F‘ 

Gj^CO) = 0 , G^C®) = 

0-^2 ■ C26+Xi)£’G2 + Cl+Xi)Fg»-26 F‘g2 = -2^’F" + % 

G 2 CO) = U = G2C«) 

where X^, d^ and are given in table 1. In 

the above equations the superscript i on F, and Gg have been 
omitted for convinience of writing. 


(3.35a) 

(3.35b) 

(5.36a) 

(5.36b) 


5.4 Skin Friction and Displacement Thickness 

The skin friction in non-dimensional form may be written 
as 
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_T p-1/2 3U 


P ut 


3N 


N=0 


Making use of (2.11a), ( 5 . 5 ) and ( 3 .I 5 ) we have for skin friction 


^ r^' u| 

— = [fi.(O) ^ r-1/2 ^ 3 + ..](3.37) 

P m=£,t,d,v ° 

To obtain the expressions for the first order and second 
order displacement thickness, to start with the basic definition 
of the displacement thickness 6 * defined as 


/ Ur^dii=J°ux^dn 

where n^ is taken as the outer limit of inner boundary layer. 
The inner expansion (5»58) may be written as 

/ ^’Cs,R"^/%)rjdN=/"[UCs,R*^/%)-uCs,N)]rj dN 

0 0 ° 


(5.58) 


( 5 . 39 ) 


The displacement thickness itself is e 2 g)anded as 
6*Cs,R“^/ 2^ = R-^/2 + p-1 + ... 

Substi luting outer expansions for U(s,E%) and inner expansions 
for n(s,E) and , we hasre 


and 


r [1 

0 


u^Cs,N) 
Ui Cs, 0 ) 


■3 dN 


( 5 . 40 ) 



’^2(s>0) _ Uj(s,N) 

UlCs,0) I ^u7Cs,0) 


U2Cs,N) - Nuj^^Cs,03 
U 2 CS.O) 


■3 dN 


2Ui(s,03 ^1 


5^ 

[ 1 - f N {1 - 
0 


Uj^Cs,N) 

UiCs, 0 ) 


-} dN3 





is 


Substituting the siinilarity -variables, the expression for the 
displacement thickness is 


5 * = 


R-V2 + R-1 ^ 


.Cd) 


r3 U, 


r5 Ui 


25 B! ^ j cos 6 

+ — L 6^ + K 5^^-* + 




where 


(3.42) 


6 = 


lim Cn-f) 

n ^ » 


(5.45a) 


W- lim 

T\ -> » 

,1 2 „(£), .. 6 
( j n - F*- •' Cn3)+^ 

(3.45b) 

= lim 

n -»• » 

(i^2 ^ 

( 3 . 430 ) 

(<^5= lim 

[n - Cn)3 - 5 

(5.43<i) 

n -»■ “ 



M . ii„ 

2 

[i^2 ^ pCv)^^33 „ i_ 

(5.45e) 


n -»• 

5*5 Heat Transfer and Eecovery Factor 

The -total heat transfer rate at the wall in non-dim^sional 

form 

q 


= .„-l r-1/2 W 

3N 


N=0 


P Uc 

Using equations (2.10), (3*4) and (5."l4) 's't® gat 

3 G' (0) + ...]/>^ 

= .a“^ P UjCs.O)Ct„-HpR"^/^[qi + R**^''^ q2+... 


0 c T U 
^ p c c 


]/m (5.44) 
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where = gi(0) + E g’(0) (545) 

q2 = I + E G^'^^’CO)] (346) 

nF£,t,d,v,e 


The temperature of the wall for whidh. there is no heat 
transfer at the surface is known as recovery temperature t^ . 
Prom equation [3-40] with q = 0» and from the expansion for 
recovery temperature 


T 


+ R 


1 

"^t 


where t and t are first and second order recovery temperature, 

we hasre expressions for t and t as 

^2 


tj, - Hi 


u: 


gi(0) 


(343) 


and 


S. . [g‘CO)GW’co).g;co)G^’"^'co)] - 

= I K — = ::: — — „t ' frn r~ 

u\ EF^,t,d,v g^^CO) ^1*“ ^ 

where i = - H‘ »^/U? 

© 2^ A 


(544) 


(545) 


Recovery temperature is generally expressed in terms of dimen- 
sionless recovery factor defined as 


r^ » 2(t^ - H^)/a^ 


1 

r^ + R~^ r^ + ... 


(546) 


where r^ and r^ axe the first and second order recovery factors 







and axe given "by 


= - 2g^C0)/g;>C0) 


(5.47) 


^2- I , 

m=£,t,d,v 


{2 Bjg^C0)Gj“^’C0) - g’CO) C0)3/g^2(03} 


+ 2 Bg 


C3.48) 







CHAPTER 4 


RESULTS 4KD DlSCUSSIOfS 


The equations for the first and second order houndarj layer 
problems (5*30 - 31 ) and (3.34 - 36) have been integrated numerically 
by Eunge-Hatta-Gill method on I.B.M. 7044 at Eanpar. The equations 
have been solved for vaxious values of suction parameter C in the range 
of -0.6 ^0.6. The ranges for other parameters ■are -7*5 ^ 

A^ <_ 4» B = 0, 0.5, and 1 and a = 0.7, 5 and 5. Resalos for second 
02?der contributions to skin friction P"(d), heat transfer G^(0) and 
recovery factor T2 are displayed graphically. The results for G^(0) 
can be obtained from those for G^(o) and 

For longitudinal curvature problem, solutions are sho'ftn. in 
Figs. (2a-c) and (3a-c). Fig, (2a) shows skin friction f"(0) verses A^ 
for C s= 0.5, and -O.5, g = O.5 and a = 0.7. For < 0, we observe 
many singularities. The location of first singularity for C = 0 (without 
suction or injection) is also shown in the same graph by a vertica,! 

3XT01! (Afzal and Oberai 1972). The first critical value (approximately) 
of for C = 0 is at A^ = ~ 3*1, for C = 0.5 (suction) at -3-5 ond 
for C = - .5 (injection) at -2.72. The effect of suction on locations 
of singularities is to make the critical values of parameters An more 
negative when compared to no suction case, i.e., the singularity is 
shifted towards left. For the case of injection the effect is otherwise. 
It may be noted that the effect of suction on locations of singularities 
is, qualitatively similar to those of favourable pressure gradients 
(Raisinghani, 1975)- i.e., for a given suction if g increases, the 





critical values move towards left. Pigs. (21) and (2c) show the effects 
of suction and injection on heat transfer G-]^(g) and recovery factor r^ 
respectively. Here the approximate location of first singularity for 
G = 0 is at = - 2.75> G = 0.5 at = - J.O aid forC = -- 0.5 
at A^ = - 2 . 32 . It is to he observed that for a fixed value of suction 
parameter C, the singularity in G^(0) and occurs s.t lo'/er valuas of 
negative ^ 'fJh.en ccmpared tothe corresponding skin friction results. This 
is due to the fact that the eigen value for the homogeneous energy 
problem is lower than the corresponding momentum problem. 

It has been pointed out by Raisinghani (l975) that the second 
order contributions for A^ < 0 sxe not of much practical importance and 
therefore detailed results have been computed for A^ > 0, and are shown 

in Pigs. (3ajb,c). Pig. (3a) shows the effect of siiction and injection 
on second order skin friction and displacement thickness. It is seen 
that for A^ =0 and for a given 3 , suction increases (in magnitude) 
the skin friction vfhile injection decreases (in magnitude) it. Pig. (3b) 
shows that for o = the second order heat transfer G^(0) decreases 

with suction and increases with injection. However for high Pranitle 

numbers (c? = 3> 5)j the trend of heat transfer cuive is reversed after 
a certain value of injection i.e., it starts decreasing. 'The rate of 
decrease beccmes more pronounced as a increases. Pig. (3c) shows results 
for recovery factor As C increases, the magnitude of recovery factor 

decreases monotonically. 

For transverse curvature problem, the results are shoTO in 
Pigs. (4a,b,c) and Pigs. (5a,b,c). Pigs. (4a,b, c) show the effects 
of suction and injection on singularities in the solutions of second order 





c_<-r 

s3d.n friction, heat transfer and recovery factor respectively. Here 
the said effects a.re similar to those discuss;! in the case of longitudinal 
curvature problem and hence need no additional comments. Pig. (52') shons 
results for second order skin friction and displacement thickness. For 
B = 0, suction reduces the skin friction while injection increases it. 

For g = .5 and 1, skin friction decreases vdth aiction (at a loner ra.te 
as compared to that for g = o) while with injection it increases, 
reaches a maximam and then decieases. The nature of ttse second order 
heat transfer curves, as shown in Fig, (5b), is similar to those for 
longitudinal curvature problem. Pig. (5c) shows the effect of suction 
and injection on second order recovery factor. For g = 0, a = 5 and 5» 
the magnitude of first increases, reaches a maximum and then decreases 
very rapidly as C decreases. 

For displacement speed problem, the results are shown in 
Pigs. (6a,b5c) and (7a,b,c). Pigs. {6a., ^, 0 ) show the effect of suction 
and injection on locations of singularities for skin friction, heat 
transfer and recovery factor respectively. The behaviour of the 
singularities is similatr to those discussed earlier in the case of 
longitudinal and transverse curvature problems and hence no further 
comments are needed. 

Pig. (7a) shows variation of second order skin friction and 
displacement thickness verses suction parameter C. Skin friction 
increases monotonically with C. For a given value of C, a(= .?) and 
S(= .5)» ■fcbe skin friction increases a.s increases. The magnitude 
of second order hea,t transfer, as shown in Fig. (7b), decreases with- 
increasing values of suction and injection. For a fixed value of C, 
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B (= ‘5) i a (= . 7 ) j the heat transfer increases ^ith . Fig. (7c) sho.'/s 
results for second order recovery factor. The recovery factor increases 
inonotonice.lly as C decreases. For B = 0, o = 3 and increases, 

reaches a maxiimm and then decreases very rapidly as G decreases. For a 
given Cj 0 (= . 7 ) and &(= , 5 ), the recovery factor increases virith A^. 

For external vorticity problem, the results are shovm in Figs. 
(8a,b,c). The second order skin friction, as shown in Fig. (Sa) decreases 
monotonically with suction parameter C. For a given value of C, the skin 
friction increases os B decreases. Fig. (8b) shows results for second 
order heat transfer. For 0 = .7, the magnitude of heat transfer decreases 
monotonically with C. For 0=3 and 5, the magnitude of heat transfer 
increases, reaches a maximum and then decreases. The second older recovery 
factor Tg shown in Fig.(8c) decreases monotonically as suction parameter 
C increases. For B= 0, 0=5 the curve shows a maxima and decreases as 
injection (-C) increases. 

For temperature gradient problem, the results are shown in Figs. (9^0 
and (Sh). The second order heat transfer increases monotonically with 
parameter C. For a given value of C and B, the recovery factor increases 
wltho. From Fig. (9b), it is observed that the magnitude of second order 
recovery factor increases monotonically with C. Further for a. given v.-.,lae 
of C, 6 , the magnitude of recovery factor increases as o decreases. 
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